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This paper presents an experimental study of the fretting crack nucleation threshold, expressed 
in terms of loading conditions, with a cylinder/plane contact. The studied material is a damage 
tolerant aluminium alloy widely used in the aerospace application. Since in industrial problems, the 
surface quality is often variable, the impact of a unidirectional roughness is investigated via varying 
the roughness of the counter body in the fretting experiments. As expected, experimental results 
show a large effect of the contact roughness on the crack nucleation conditions. Rationalisation of 
the crack nucleation boundary independently of the studied roughnesses was successfully obtained by 
introducing the concept of effective contact area. This does show that the fretting crack nucleation 
of the studied material can be efficiently described by the local effective loadings inside the contact. 
Analytical prediction of the crack nucleation is presented with the Smith- Watson- Topper (SWT) 
parameter and size effect is also studied and discussed. 
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I. INTRODUCTION 

Weight saving is still a topical issue in aerospace ap- 
plications. Because of their low densities, aluminium al- 
loys are widely used in air plane structures, in spite of 
quite bad fatigue and tribological properties. As noted by 
many authors [J [2]) a very important example of applica- 
tion is riveted lap joints structures, where contact fretting 
stress, due to vibrations during the flight, can initiate 
cracks at the hole edge which may propagate due to the 
overall structure fatigue stress. In the worst case, propa- 
gation will lead to failure. Authors 0, Q also report par- 
tial slip in rivet structures. This shows the importance 
of investigating fretting crack initiation boundary in par- 
tial slip conditions of concerned material. In addition, 
in industrial application, where variable surface qualities 
can be found, one major issue is the impact of surface 
roughness on the crack nucleation process. Coupling ex- 
periments and modelling this paper focuses on this topic 
in order to predict the cracking risk on a 2024T351 alu- 
minium alloy under fretting loading and taking into ac- 
count the roughness impact. To evaluate the roughness 
effect on a simple fretting test, a cylinder/plane configu- 
ration with three different fretting pads in terms of rough- 
ness, was used. After a brief description of fretting test, 
material studied and surface roughness morphologies, the 
crack nucleation boundary at S-IC* fretting cycles under 
partial slip condition is defined for the smooth contact 
case. Two unidirectional surface roughnesses are succes- 
sively analysed and the corresponding crack nucleation 



boundary is determined. Introducing the concept of ef- 
fective contact area, the different experimental results 
are rationalised through the introduction of a crack nu- 
cleation master curve. Finally, applying the SWT multi- 
axial fatigue criterion, prediction of the crack nucleation 
is intended taking into account the size effect. 



II. EXPERIMENTAL WORK 

Fretting is defined as a small oscillatory movement be- 
tween two bodies in contact which induce a relative dis- 
placement between the two surfaces. Considered as a 
plague for modern industries, it is a very complex prob- 
lem involving many aspects as tribology, contact fatigue 
mechanics and material science. To reproduce such phe- 
nomena, different fretting test devices have been devel- 
oped 0,11. 



A. Material and specimens description 

2024T351 is an aluminium copper magnesium alloy 
which has been solution heat treated, control stretched 
and naturally aged. Chemical composition and mechan- 
ical properties are listed in tables U and |TT1 respectively. 
R-Po.2% denotes the classical yield stress defined at 0.2% 
of strain ; E the young modulus, u denotes the pois- 
son's ration and ad the fatigue limit. Specimens were 
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TABLE I: Chemical composition of Al 2024T351 (weight %) 
%Cu %Mg %Mn %Fe %Si %Cr %Ti 
4.4 1.45 0.62 0.20 0.15 0.01 0.03 
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TABLE II: A summary of the mechanical properties of the 
2024T351 alloy used for fretting tests. 

Rpo.2%(MPa) E(GPa) v crd(MPa) 



325 



72.4 0.33 



140 



extracted from a laminar plate of 25 mm thick sup- 
plied by Pechiney Industries. After machining, the sur- 
face in contact were carefully polished until mirror state 
(Ra around O.OS/xto). Displacement direction of the fret- 
ting tests is parallel to the L direction of the material 
along which grains are elongated(see fig. [l]). Grain size 
in the T direction was determined to be 150/im by Elec- 
tron Back Scattered Diffraction (EBSD) mapping, using 
a line intercept measurement method. 
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FIG. 2; Detail of fretting test device. 




FIG. 1: Al 2024T351 microstructure of the sample surface. 
Optical micrograph, after a Keller reactive attack, show- 
ing elongated grain structure in the rolling direction. Inter- 
metallics particles appear in dark. 



is made by cutting the specimen in the middle of the 
contact zone and polishing this surface. Crack initiation 
is investigated through optical micrography and images 
are recorded if any damage is visible. The crack angle 
with respect to the surface, and then crack depth below 
the surface are measured. Such a cross section of a fret- 
ting crack observed by optical micrography is presented 
in figure [3l 




B. Fretting experiment 

The experimental setup used in this study is based on 
a fretting device rigidly mounted to a servo hydraulic 
test machine (see figure [2]). The counter body is made of 
7075T6 aluminium alloy ; the cylinder radius is 49 mm 
and the cylinder length L is 4.4 mm. In cylinder plane 
configuration, it is helpful to define normalised loads P 
and Q with respect to the contact length (eq. [T|) , which is 
perpendicular to the sliding direction. The linear normal 
load range is 100 to 1000 N/mm, equivalent to a nominal 
maximum hertzian pressure from 220 MPa to 700 MPa. 



and g = ^ 

Li 1j 



(1) 



Once specimen is fixed and carefully aligned with the 
counter body, fretting solicitation is applied during 5.10^ 
cycles by monitoring the normal load imposed to the sur- 
faces (P) and the displacement amplitude {6) or the tan- 
gential load amplitude induced by the contact {Q). For 
the oscillatory parameters the stared versions {5*^Q*) 
denote their amplitude. After the test, an optical micro- 
graph of the fretting scar is recorded and a cross section 




'center of the contact 



FIG. 3: Optical micrograph of the cross section below the 
surface. 



C. Surface quality 

Relatively few literature can be found on the effect of 
roughness in fretting damage. Although a detrimental 
effect of the roughness is expected, some authors have 
reported a great benefit in term of lifetime can be ob- 
tain by a special laser surface texturing In order to 
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study the roughness effect on the fretting crack nucle- 
ation conditions, three different cyhnders named Rl, R2, 
R3 with increasing roughness, are used. Each counter 
body is obtained by lathe machining so that roughness 
is assumed unidimensionaL R3 is machined with rela- 
tively high speed advance of the tool, R2 is obtained 
with a lower one and Rl is polished on the machine af- 
ter as fine as possible machining. Figure [3] presents the 
three profiles of the fretting pads and one of a specimen 
for comparison. Respective roughnesses are measured on 
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FIG. 4: Comparison of the 3 fretting pads roughness profiles. 



a tactile profilometer and the corresponding results are 
gathered in table IIIII 



TABLE III: Mean and maximal roughness of cylinders and 
specimens 

Rl R2 R3 Samples 

Ra(/^m) 0.11 0.60 0.75 0.04-0.08 
Rtium) 1.05 3.10 3.15 0.082 



The main fretting damages can be rationalised in a Ma- 
terial Response Fretting Map (MRFM) involving non- 
degradation, cracking, particle detachment in a map with 
normal load versus relative displacement. The fretting 
description proposed by Vincent et al. Q is based on 
these two sets of fretting maps (see figure [5]) • When the 
relative displacement of the surface is small, partial slip 
prevail and induce crack initiation, although at higher 
values of the displacement, gross slip settles and impose 
wear of the surface. Since in our case, partial slip is 
assumed, very small displacement amplitudes have been 
maintained to study the crack initiation. 

running condition fretting map (RCFM) 




material response fretting map (MRFM) 




nt ampHtude [/im]-^^ 



FIG. 5: Running condition fretting map which combines the 
fretting regime analysis (RCFM) with the material response 
(MRFM). 



B. Experimental fretting regimes determination 



III. TRIBOLOGY ANALYSIS 

A. Fretting regime analysis and material response 
maps 

Two essential sliding conditions can be defined: Par- 
tial slip regime (PSR) characterised by a closed elliptical 
fretting loop, with a composite contact of a sliding and 
a sticking zone ; and Gross slip regime (GSR) identified 
by a quadratic dissipative fretting loop, with full sliding 
occurring over the entire contact zone The transi- 
tion between one regime to another defines the Mixed 
fretting regime (MFR). Identification of this regimes in 
carried out by plotting a Running Condition Fretting 
Map (RCFM). Waterhouse [8] first highlighted a corre- 
lation between the sliding regime and fretting damage. 



To determine the fretting regimes, variable displace- 
ment method [l3|, has been applied for different val- 
ues of the normal load. In each test, the normal load 
is kept constant. The relative displacement start to a 
very low value, imposing partial slip (with Q* < HtP)- 
When stabilised conditions are reached, d* is increased 
and maintained constant until new stabilisation is reach. 
6* is increased this way, step by step, until the con- 
tact has overcome the sliding transition conditions (i.e. 
Q* — iitP)- Precise value of the transition is then cal- 
culated by computing the value of the energy sliding cri- 
terion A = Ed/Et ; with Ed the dissipated e nergy o f 
the corresponding cycle and Et the total energy [llf, [13] ■ 
In the cylinder/plane configuration, it has been shown 
that the sliding transition is associated to a discontinu- 
ity of ratio A (see figure [6]). Figure [7] present the slid- 
ing transition obtained with the variable displacement 
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FIG. 6: Illustration of variable displacement method (see the 
text for details), Fn=750N. The ratio Q*/P (left scale) in- 
creases until stabilised gross slip condition is reach. The en- 
ergy ratio A is computed and exhibits a discontinuity at the 
transition. 



method, for a cylinder/plane contact of A12024T351 ver- 
sus A17075T6. On this figure, results are presented ver- 
sus the relative displacement amplitude d* . However, it 
does not appear as the more convenient representation 
when different series of results are compared. Indeed, 
the value of 6* depends on the fretting setup compliance 
[T^ . Since we will use 3 different fretting pads and dif- 
ferent contact compliances, the displacement amplitude 
is adequately substituted by the linear tangential force 
amplitude Q*(N/mm). Note that the tangential force is 
not affected by the apparatus compliance. In the next, 
except if specified, results will be presented using the 
tangential force. 
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FIG. 7: Experimental determination of the sliding transition, 
cylinder/plane contact, 2024T351 vs 7075T6, R=49mm. 



C. Coefflcient of friction in partial slip 

Under partial slip conditions, slip takes place in two 
symmetrical regions a > \x\ > c which surround a central 
stick zone < c. In this central zone, from Coulomb's 
law we have q{x) < x p{x). Because of the stick zone, 
the friction coefficient obtained at the sliding transition 
cannot not used. Instead, the friction coefficient in par- 
tial slip regime /ips can be determined by the investiga- 
tion of the scars, after the fretting tests. From Mindlin 
cited by Jonhson [l^ we have: 



c = a \ 1 



Q* 



1/2 



From which: 



MPS 



Q* 



1 



(2) 



(3) 



From examination of fretting scars, the contact width 
2a and the stick zone width 2c (defined on figure [8]) are 
measured for a large range of displacement values and 
the friction coefhcient in partial slip regime is calculated 
in each case. This was done for a constant normal force 
Fn = 1400A'' (P = 3187V/mm) and result is assumed to 
be true wathever the value of P is for the investigated 
experimental conditions. Figure [5] summarise the results 
obtained by this methodology. The evolution toward 
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FIG. 8: Experimental scar after a fretting test, Fn=1400N 
and Ft=500N. Definition of stick and slip zones width. 



gross slip as displacement (and so tangential force) in- 
crease is well observed. Calculations of the friction co- 
efficient in partial slip regime with equation ([3]) show 
clearly a constant evolution, very close to the value at 
the transition. Then we can postulate that the coeffi- 
cient of friction under studied partial slip conditions, is 
identical to the value defined at the transition. 



(4) 



IV. 



CRACKING ANALYSIS IN PARTIAL SLIP 
CONDITIONS 



The experimental objective in this part is to determine 
as finely as possible for each studied normal force, the 
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FIG. 9: Determination of the friction law in partial slip 
regime, Fn = 1400Af. 



critical tangential force associated to the crack initiation 
after 5.10"* cycles under stabilised partial slip conditions. 



A. Methodology 

For each normal load, kept constant during a test, sev- 
eral fretting tests were carried out as detailed in IIIBl 
maintaining a constant tangential force (Q*) by mon- 
itoring the displacement amplitude (6*). The critical 
value (Q*) is precisely determined thanks to cross sec- 
tion observations of specimens cut after fretting solici- 
tation. Cracking is assessed by optical inspection, the 
smallest observable crack being around 5fim long. To 
bracket the crack nucleation threshold, dichotomy is ap- 
plied to select the values of Q* About ten test are 
required to reach a 5N accuracy in the determination of 
Q*. Measurement of the crack lengths allows to extrap- 
olate Q* even if dispersion can sometimes disturb the 
method. This methodology is illustrated in figure [TUl 



B. Results with a smooth contact 

The described methodology is applied with Rl counter 
body for 5 normal forces and crack nucleation bound- 
ary is determined (figure [TT|). In addition. Several other 
conclusions can be derived from experimental observa- 
tions. First of all, the crack nucleation is always located 
near the trailing edge of the contact, which is in agree- 
ment with former fretting studies in partial slip regimes. 
When crack initiation is detected, the crack angle is mea- 
sured and is found to be lower than 35 degrees for all the 
loading conditions investigated. Referring to the work of 
Lamacq [3] and Dubourg and Lamaca [iq , cracks seems 
to growth under mode II by maximal shear amplitude. 
Figure [11] shows that the threshold tangential load is 
found independent of the normal force, giving: 
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FIG. 10: Schematics of the crack nucleation threshold deter- 
mination: a) fretting test conditions ; b) results in terms of 
crack length. 
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FIG. 11: Experimental fretting crack nucleation boundary, 
cylinder/plane smooth contact (R=49mm). 



V. ROUGHNESS IMPACT 

In this section, results with rough surfaces are pre- 
sented and compared. As expected, significant gaps are 
observed between each crack nucleation boundary. The 
effective contact area concept is introduced in order to 
account for the roughness. Finally, a new representa- 
tion is introduced which rationalises the different crack 
boundaries obtained for different roughnesses. 
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A. Results with rough contacts 

Mean and total values of fretting pads roughness were 
presented in table lTTTl The same methods as used in para- 
graph IIII Bl for determining the sliding transition and 
in paragraph IIVBI for establishing the crack nucleation 
boundary, were used with the fretting pads of roughness 
R2 and R3. The corresponding results are shown in fig- 
ures [12] and [121 On fig. [T^ it can be seen that the sliding 
transition does not depends on the roughness value. This 
could be expected since gross slip condition is imposed at 
the transition, which induce wear in the contact. Rough- 
ness is then blurred in a few cycles; and the contact in- 
terface tends to a unique structure. On the contrary, the 
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FIG. 12: Comparison of sliding transitions obtained with 
smooth and rough cylinder. 



crack initiation, under stabilised partial slip condition, 
strongly depends on the surface roughness (figure [T^ . 
As expected, a higher roughness leads to a lower tan- 
gential threshold load. Fig [T3| also shows a slope in the 
experimental boundaries obtained with R2 and R3. This 
influence of the normal force on the tangential thresh- 
old load was not observed with Rl. Those differences in 
tangential load threshold can be explained in considering 
intrinsic contact parameters. Indeed, the observation of 
fretting scars f figure [H]) . clearly shows that for identical 
loading conditions, the fretting pad surface in contact 
varies with the surface roughness. This implies different 
loading conditions inside the contact, depending on the 
area effectively involved in the contact. 



B. Definition of the effective contact area 

First developed by Greenwood [l6|, the effective con- 
tact area S'e// can be simplified in our case. Considering 
the polished surfaces of the specimens and the unidimen- 
sional fretting pad roughness, Sej / is easily defined as the 
sum of the micro-rectangles of the pad indentation. This 
definition can be one more time simplified by 2a x ie// 
where Le/ / is the sum of the micro- length of indentation 
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FIG. 13: Experimental determination of crack nucleation 
boundaries with three different roughness. 



on the medium line of the contact f figure [T5)) . This rather 
simplistic definition is only valid for well aligned contact. 
Next step is to quantify the influence of -F/v and Ft on 
Le/ / ■ Effective contact lengths were measured from opti- 
cal images of fretting scars after the tests. Results clearly 
show a negligible influence of tangential load so that Le/ / 
appears mainly monitored by the normal load. A linear 
dependence of S'e/ / with Fj^ is observed ; the quasi hori- 
zontal curve obtained for Rl shows that S'e/ / = in the 
investigated loading range (figure [16]). 



C. Unification of the fretting crack nucleation 
boundary 

Two new loading parameters Pgy / and Q^j f , respec- 
tively the effective normal load and the effective tangen- 
tial load are introduced : 

Fj\[ Ft 
Peff = 7 and Qeff = (6) 

Fig. [16] gives the evolution of Le/ / with Fn ■ Assuming 
a negligible impact of the tangential force on Le/ f , the 
values of Peff and Qeff can be computed. All the fret- 
ting tests defining the three nucleation boundaries can 
then be plotted in terms of the intrinsic contact loading 
parameters Pef / and Qef / . The corresponding represen- 
tation is presented on figure [TTl On this figure, compared 
to the fig. 1131 the crack nucleation boundary under par- 
tial slip conditions appears independant of the roughness 
value (for the experimental conditions investigated). Ac- 
cessing intrinsic contact loading parameter was successful 
in correlating the three experimental boundaries. One 
point further reinforce this correlation: in fig. 1131 the 
slope with R2 is larger than with R3 even if the rough- 
ness is lower. Figure [T6l shows that the slope in the ie// 
evolution with Fm is larger for R2 than for R3. This re- 
sult might seem incoherent with the fact that zero slope 
is expected for a perfectly smooth contact. However it 
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FIG. 14: Examples of optical micrography of specimen sur- 
faces after a fretting test, a) pad Rl, Fjv = 1400A'^ and 
Ft = 500iV b) pad R2, Fn = SOOiV and Ft = 400iV c) 
pad R3, Fjv = 600iV and Ft = 450iV. 
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FIG. 15: Definition of the effective contact area 5e// (see the 
text for details). 
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FIG. 16: Experimental determination of Le// evolution with 
normal force (pad length : L — AAmm). 



of P observed with the two rough contacts on fig. [121 
This effect may thus be related to the material response 
versus the indentation of the asperities of the counter 
body. 



could be explained by a different response of roughness 
to indentation induced by the quite different shapes of 
the roughness profiles shown on fig. |4l This outlines 
the difficulty to rationalise a surface morphology by only 
considering the two parameters Ra and Rf. In any case, 
the larger slope observed for the L^f f-R2{FN) curve per- 
fectly accounts for the good correlation observed on figure 
[T7]and highlights the physical reliability of the approach. 
To summarise fig. I17[ it is shown here that the crack ini- 
tiation of the 2024 alloy studied depends on the real lo- 
cal stress state {Peff, Qeff), which cannot be calculated 
without considering the real contact length L^ff, in case 
of a rough contact. It is also shown here that considering 
the effective loading parameters have removed the effect 



VI. FRETTING CRACK INITIATION 
PREDICTION WITH SWT PARAMETER 

To predict the crack nucleation risk and formalise the 
previous experiments, the SWT multiaxial fatigue cri- 
terion is implemented and combined with a Mindlin de- 
scription of the cylinder/plane contact configuration. Ex- 
plicit formulation will provide a physical justification to 
a non dependence of the crack initiation boundary with 
respect to the normal load. Size effect is also introduced 
and implemented. It should be noted that alternative ap- 
proaches like combining classical fracture mechanics and 
notch fatigue can also be used to tackle this problem. 
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FIG. 17: Experimental fretting crack nucleation boundaries, 
with various pad roughness, plotted in Pef f /Qef f represen- 
tation. 



Such an approach was recently carried out by Ciavarella 



A. Theoretical background 

The SWT critical plane approach was first used to 
tackle fretting problems by Szolwinski and Farris [if. 
They used a fatigue approach developed by Smith et al. 

which was developed to account for a mean stress or 
strain effect in classical fatigue. According to the Szol- 
winski and Farris approach, the crack initiation in fret- 
ting occurs on the plane where the product of the nor- 
mal strain amplitude ta , and the maximum normal stress 
Cmaa; IS maximum. The SWT parameter F can then be 
expressed as: 



X ea = 



E 



-i2N) 



2b' 



'J2Nf 



(7) 



where a'^ is the fatigue strength coefficient, 6' is the fa- 
tigue strain exponent, e'j is the fatigue ductility coef- 
ficient, c' is the fatigue ductility exponent and N the 
number of reversal to failure. These constants are also 
given by Szolwinski et al. 19] and listed in table HVl 
To normalise the crack nucleation risk, the following 



TABLE IV: SWT parameters for Al 2024T351 from 19] 
o-f(MPa) b' e'f c 
714 -0.078 0.166 -0.538 

scalar parameter dswT can be introduced : 



If dswT is greater than or equal to 1, cracking is likely 
to occur. If dswT is less than 1, there is no risk of crack- 
ing. By expressing analytically Umax and Ca, a literal 
expression of the cracking risk dswT can be obtained. 
Fouvry "12] and Fridrici 20] showed that the SWT multi- 
axial fatigue criteria application on fretting contact gives 
a maximum crack initiation risk a the trailing edge of the 
contact. This is in agreement with our experimental ob- 
servations. We can therefore express the surface stresses 
at the contact border in partial slip with Mindlin theory. 
Assuming an elastic plane strain 2D state, we use the 
general equation [2l| : 




Wo 



i V g ' a J 

fJ-Po 



(9) 



Where af^ and ajj denotes stresses in a cylinder plane 
contact submitted to normal load and tangential load 
respectively. Note that the subscript ij can denote any 
spatial direction. At the trailing edge of the contact (at 
X = a and y = 0) we have : 
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<(i,o) 
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(1,0) = 
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and 
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' xy 



2Aipof -2A*Pov/(f)^-l 


^^xa*,(f,0) 
(f,0) = 



(10) 



which gives : 



(Jxx{a,Q) = 2fipo^ 1 - (^^^ 



(11) 



This equation may be written in a more convenient form 
by combining it with ^ to give : 



o-xx{a,0) = 2/ij3oA 



In cylinder /plane configuration we have po = ("9^) 

and experimental analysis fsee llll C|) showed that ^ps 
/i = 1.1. Combining these elements we get : 



(12) 

1/2 



^^^^(a, 0) = 2 



IJ.E*Q* 
ttR 



(13) 



The stress state at the trailing edge of the contact can be 
considered as quasi uniaxial. From Hook's law we deduce 
at point (a,0) : 
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On the loading path of a fretting cycle, we have at point 
(a,0) : 



AE 



(15) 



The SWT parameter is finally simplified to the following 
expression : 



r = 



E 



(16) 



We consider the critical crack initiation limit [dsWT — 
1) so that the equivalent critical crack tensile stress, with 
respect to SWT criterion, is expressed as : 



1 - 2j/2 - J/3 



\b'+c' 



(17) 



Finally from equation (|13p we determine a critical linear 
tangential force amplitude associated to the crack initia- 
tion condition : 
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ty 

averaged loading path 

transposition to multiaxial fatigue criteria 

FIG. 18: Schematic representation of the radial shape process 
volume averaging. 



vri? {a'jn2Nr + ^',e'j{2Nf 

- liiE; 1-2^/2-^3 (18) 

The critical linear tangential force amplitude is propor- 
tional to the cylinder radius, a function of the elastic and 
fatigue properties of the material and inversely propor- 
tional to the friction coefficient. It must be pointed out 
that Q* does not depend on the normal force which ac- 
count for the observed negligible impact of the pressure 
field on crack nucleation. Performing a numerical appli- 
cation leads to Q* — 106iV/mm, value to be compared 
with the experimental threshold of 240 N/mm. This dif- 
ference could be related to the large stress/strain gradi- 
ents which exists in and below the contact region. 
Although this approach gives a theoretical justification of 
the normal load non dependance, it cannot predict the 
exact threshold tangential load. Indeed, fretting crack 
initiation behaviour may also depend on the severity of 
the stress/strain gradients, and not only on the maximum 
values. First introduced by Fouvry et al. (T^. [2^. the 
crack nucleation process volume was successively adapted 
by Araujo and Nowell [1^ and more recently used by 
Naboulsi and Mall [l^]. In this approach, the compu- 
tation of a representative SWT critical plane parame- 
ter must be conducted from an averaged stress strain 
state, defined over an intrinsic processing material vol- 
ume. Hence, the stress gradient effect is indirectly con- 
trolled through the averaging process. 

B. Crack nucleation process volume approach 

To compute the SWT parameter using the process 
volume approach, a cross section is discretised and par- 
tial slip stresses and strains corresponding to the experi- 
mental conditions, are computed from Mindlin solutions. 



Process volume averaging is then applied considering a 
radial shape volume as shown on figure [181 The SWT 
crack initiation parameter is then computed from these 
averaged values which gives a mean SWT value on the 
process volume. As discussed in [l^], SWT could be com- 
puted on each point, the averaging procedure occurring 
thereafter ; but the result would not be very different 
and the first method, which is faster, has more physi- 
cal meaning. From these computations, the critical lo- 
cations where cracking risk is maximum, is located at 
the contact surface in the vicinity of the trailing edge. 
This matches the experimental observations seen in sec- 
tion IIVI However if the crack location at the trailing 
edge is well predicted, the SWT parameter does not pre- 
dict the observed crack initiation angle (<35 degrees). 
All the values are closed to 90 degrees, perpendicular to 
the maximum principal stress. Indeed, mainly dependent 
on the tensile stress and strain state, the SWT param- 
eter is unable to describe the shear mechanisms which 
presently control the crack nucleation process. The pro- 
cess volume has been identified for a linear effective force 
of 380 N/mm. Figure [TOl displavs the surface distribution 
of the cracking risk dswT defined in equation ([8]), for dif- 
ferent averaging length scale values. Experimental cor- 
relation with the crack nucleation boundary is obtained 
for a r = SOfim process volume radius. This value match 
very well the mean grain size radius of Ibfim, and con- 
sequently confirms a correlation between the mechanical 
crack nucleation process volume and the microstructure. 
With the process volume radius fixed, by varying normal 
load in the computation analysis, one can draw the ana- 
lytic crack initiation boundary predicted with SWT pa- 
rameter. This is compared to the experimental boundary 
on fig. [20] and plot in the effective loads representation. 
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FIG. 19: Application of the averaging method in calculation 
of the cracking risk at surface with the SWT parameter. 
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FIG. 20: Comparison between experimental and predicted 
crack nucleation boundary with a 80/im size effect. Fret- 
ting contact: A12024T351 plane versus A17075T6 cylinder pad 
(R=49mm). 



C. discussion 

Obviously, figure [20l confirms that SWT computations 
with a.r — SOfim process volume radius allows to predict 
the studied 240 N/mm linear normal loading. Although 
a weak influence of the normal load is observed, it is 
nevertheless deduced that the application of a singular 
process volume radius is ineffective to quantify the pres- 
sure impact. Indeed, in the SWT formulation, pressure 
effect have tendency to reduce the cracking risk which 
is based on tensile stress, although the contrary is ob- 
served experimentally. Thus the physical meaning of the 
SWT parameter to predict the 2024T351 fretting crack 
initiation appears discutable. Other criteria, more repre- 
sentative of the physical behaviour of the studied alloy, 
should be studied involving elasto-plastic analyses of the 
stress/strain loading path. 



D. Identification of a safe crack nucleation process 
volume 



The results presented above show that analytical cal- 
culations coupled with the process volume approach and 
the SWT criterion are limited for predicting the experi- 
mental crack nucleation boundary. However, due to the 
simplicity of its formulation, it is widely employed in in- 
dustrial component conception. An alternative method 
to reach a safe prediction for crack initiation in contact 
submitted to fretting loading, is presently introduced. 
First, the contact pressure range must be defined. Then 
the process volume is evaluated for both Pmin and Pmax ■ 
Assuming a monotonous evolution between Pmin and 
Pmax, the safe crack nucleation process volume Ts is de- 
fined as : 

rs = Min[r{P,nax),r[P,nin)\ (19) 

In the present case, a positive slope with the pressure 

Experimental 
crack nucleation 
boundary 

/ 

Analytical 
prediction 




Tangential load amplitude Q* 



FIG. 21: Definition of the safety domain predicted by ana- 
lytical computations of SWT parameter with the safe crack 
nucleation process volume rs. 



increase has been identified. Therefore we can directly 
assume that the safe crack nucleation process volume 
correspond to the process volume defined for Pmax, the 
maximum pressure loading. Setting Pmax to lOOON/mm 
leads to a critical value of r = 70/im which can be used to 
define the safety domain (see figure [^T|l . It is important 
to mention that the obtained process volume is very close 
to the value previously defined for the minimal pressure. 
Again, it confirms the stability of the approach and sup- 
port the idea of a reliable grain size averaging approach. 



VII. CONCLUSION 

Fretting crack initiation under partial slip conditions 
is investigated through an experimental approach and 
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the effect of rouglmess is studied quantitatively. Alu- 
minium alloys contact, under cylinder/plane configura- 
tion is studied in order to tackle crack initiation in riv- 
eted assemblies. Three fretting pads are used with tree 
different surface qualities, with a unidimensional rough- 
ness morphology. 

First, the running condition fretting map is determined 
and the friction law is quantified, the friction coefficient 
in partial slip is demonstrated to be equal to the friction 
coefficient at the transition /Xj = 1.1. The sliding tran- 
sition appears to be independent of the roughness value, 
which could be explained by a rapid blurring of the in- 
cipient surface roughness under gross slip conditions. 

Further work has been undertaken to determine the ex- 
perimental crack nucleation boundaries with respect to 
the surface roughness. Results have shown that a critical 
tangential load could be found with a smooth contact, 
although a weak pressure effect is observed for the two 
rough contacts. As expected, crack initiation appears 
very sensitive to the surface quality, a higher roughness 
leading to a lower value of the tangential force needed for 
crack initiation. A contact effective area formulation is 
introduced in order to access the intrinsic contact load- 
ing parameters. Linear normal and tangential effective 
loads Peff and Qeff are computed and a new represen- 
tation is drawn to describe the crack nucleation bound- 
aries. This is successfully applied to correlate experi- 
mental results and various crack nucleation boundaries 
are unified independently of the surface quality leading 
to Qlfj:^ = 2407V/mm. 



An analytical justification of the weak dependence of 
Q* to the normal load for a smooth contact has been 
conducted. SWT criterion has been formulated in or- 
der to predict the crack nucleation threshold. The non 
influence of the normal load on crack nucleation thresh- 
old is justified although the numerical application ap- 
pears to be on the safe side. A radial process volume ap- 
proach is then introduced to take into account the severe 
stress/strain gradient in the contact. A global correla- 
tion is found with the experimental results with a pro- 
cess volume diameter which match the mean grain radius 
of 75 /iTO and thus, this effect could be correlated to a 
stress/strain microstructure accommodation. In spite of 
this, the pressure effect on the crack nucleation bound- 
ary, is not in agreement with experimental results. To 
avoid a large underestimation of the cracking risk, a safe 
crack nucleation process volume identification approach 
is introduced. It consists in defining the smallest vol- 
ume through the investigated pressure and contact size 
ranges. 

However, some critical aspects like the pressure loading 
impact or the incipient crack orientation suggests that 
other fatigue approaches must be considered and elasto- 
plastic stress/strain analysis are required. 
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